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Abstract. 

The traditional design of trigger algorithms for GRB experiments requires the specification of the background and burst 
samples in terms of acquisition times that are of fixed duration and of fixed elapsed time from each other. One such set of 
acquisition times is required for each characteristic timescale of GRB variation that one desires to detect. One then slides 
each set through the trigger data searching for samples that maximize the signal-to-noise of the background-subtracted burst 
sample. 

Here we describe a new triggering approach in which the times at which the background and burst samples are acquired 
are allowed to vary dynamically. Two background samples bracket a burst sample. The background and burst durations and 
elapsed time between them are allowed to be free parameters, which are maximized using the downhill simplex method. This 
produces great flexibility in the timescales that are available for detecting GRBs. 



INTRODUCTION 

The search for untriggered GRBs is has been an ac- 
tive field of research at least since the public release of 
BATSE data. Most recently, Kommers et al. [1] and Stern 
et al. [2] have described searches of BATSE data directed 
at revealing GRBs that occurred without being detected 
by the BATSE on-board triggers, either for operational 
reasons or because their spectral or temporal morpholo- 
gies were poor fits to the on-board trigger criteria. 

Naturally, the same considerations apply to GRB de- 
tection by HETE. The HETE mission deploys an un- 
precedentedly varied set of trigger criteria — the FRE- 
GATE DSP trigger [3] uses four timescales and operates 
in two energy bands, while the WXM XG trigger [4, 5] is 
typically configured to apply thirty or so criteria, some on 
WXM data and others on FREGATE data. Nevertheless 
the variety of GRB morphologies, and operational con- 
siderations, can result in GRBs that are not detected in 
flight. It is important to develop a strategy to mine HETE 
survey data for such untriggered GRBs. 

Typical ground searches for untriggered bursts use 
detection methods that largely mirror on-board trigger 
algorithms [1, 2]. Background and burst samples are 
specified in terms of acquisition times that are of fixed 
duration and of fixed elapsed time from each other. Each 
GRB timescale — risetime or duration — is probed by 
a different fixed choice of these parameters. Each such 
fixed set of time windows is then swept through the 
time series being probed for transient events, searching 
for samples that maximize the signal-to-noise of the 
background-subtracted burst sample. 



This scheme has the disadvantage of being rather in- 
flexible about the the timescales that are probed. This in- 
flexibility is especially troublesome when seeking weak 
signals, for which inaptly chosen burst or background 
samples may lead to a signal dilution that prevents de- 
tection. 

In this work we describe an alternative approach that 
has been quite successful in identifying extremely weak 
events. In this approach, the background and burst sam- 
ples are treated as free parameters, which are varied us- 
ing the downhill simplex method of Nelder & Mead [6] 
to maximize the signal-to-noise ratio of the background- 
subtracted burst sample. 



IMPLEMENTATION 

The operation of the code, spiffy-t rigger, is illus- 
trated in Figure 1 . 
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FIGURE 1. Illustration of the definition of the time parame- 
ters that are varied to search for transient events. 



The trigger operates as a simplified "bracket trigger" 
[4, 5], in that the background is estimated using samples 
before and after the burst sample. The background is as- 
sumed constant, so no interpolation (linear or otherwise) 
is performed to obtain the background rate during the 
burst sample. 

This restriction is not an essential feature of the 
method, but rather merely a simplification. The two 
background intervals are restricted to remain equidistant 
from the burst sample, so as to prevent the maximization 
procedure from exploiting a monotonic increase or de- 
crease in the background rate to estimate an erroneously 
low background, by driving one of the background sam- 
ples to a region of lower background without driving the 
other to a region of higher background. 

The code operates on a time-series of integer counts. 
It advances a trigger window of fixed duration through 
the time series by steps of size . It sets up a burst 

sample interval, of duration , bracketed at an elapsed 
time T g j by two background sample intervals of duration 
Zj^ , the second of which ends at time f gn( j . It calculates 
the SNR for the burst sample, assuming a background 
rate calculated by a weighted average of the count rates 
in the two background intervals. 

The SNR is computed as follows: Assume for the 
sake of generality that the two background accumulation 
times may differ, so that we accumulate ^ counts in 
the first background during an accumulation time Tj^i , 
and n^Yi coun ts m the later background accumulation 
time • Denoting the estimated background counts 
during the burst sample by jit^ , and assuming the Gaus- 
sian approximation to the Poisson distribution, it is a 
straightforward exercise in Gaussian estimation to show 
that 
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where E 2 is the variance in the estimate jit^ . 

Denote by the counts that we accumulate during 
the burst sample. Then the net signal in the burst sample 
is s = — ■ The variance in s is the sum of 
E 2 and the variance in h^ u . Triggering is essentially 
hypothesis testing, with the null hypothesis consisting of 
the assumption that the count rate in the burst sample 
is the same as what is estimated using the background 
samples. Thus the appropriate choice for the variance of 
«k u is "model variance", that is oj 2 ^ = jU^ u . Thus the 
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This is the quantity that spiff y-trigger endeavors 
to maximize. 

The code uses the simplex method to vary the four pa- 
rameters r gn( j , , T g j , and T^ u , which are viewed 
by the simplex minimization routine as continuous pa- 
rameters. A very lax convergence criterion is imposed — 
the absolute variation of the SNR must be less than 0. 1 
across the simplex — because in triggering there is no 
point in determining the SNR to great accuracy, and be- 
cause we don't want to spend many CPU cycles chasing 
noise. 

The parameter f gn( j is constrained to be later than the 
end of the trigger window in the previous invocation. 
Consequently, the arrangement of burst and background 
samples "accordions out" backwards in time from the 
current time, without repeating choices of intervals made 
during previous iterations. 

When there is no transient event in the data, the sim- 
plex will typically not wander very far from its initial 
configuration. On the other hand if there is a transient 
event, and the initial simplex includes a vertex corre- 
sponding to a configuration in which the burst sam- 
ple even partially includes the event, the simplex will 
rapidly climb the SNR slope, dynamically adjusting its 
timescales until the event is well-bracketed. 

Since the simplex does not wander far if it doesn't 
find much at the outset, it is important to ensure that 
T skip ' s no * so ma t a short event may "fall between 
the cracks" — that is, fail to have any of its constituent 
time samples included in a burst sample probed by the 
initial simplex. It is therefore a good idea to ensure that 
at simplex initialization, t^ u > for at least one of 

the simplex vertices. This ensures that every data sample 
passes through the burst sample of at least one initial 
simplex parameter vertex. 

Constraints on the time parameters are imposed by 
making the SNR function return a large negative value 
when the constraints are violated. The previously- 
discussed constraint on the parameter t ^ is enforced 
in this way. The code also uses this parameter-constraint 
mechanism to prevent intervals from encroaching upon 
each other, to ensure that , T g j , and T^ u remain 
positive-valued, and to keep all intervals inside the 
current trigger window. 

Other useful constraints that it is good practice to en- 
force are a minimum value for T g j (so that burst and 
background samples are well-separated), a minimum du- 
ration for (so as to minimize the risk of the back- 
ground nestling into a low fluctuation), and a maximum 
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FIGURE 2. GRB detected by the trigger robot on 19 June 
2001 at 15:17:01 UTC. The SNR of this event was 6.1, the 
duration was 3.7 s. This event was also detected by BeppoSax. 

duration (so as to minimize the risk of triggering on 
very long duration trends in the background). 

DEPLOYMENT 

spif fy-trigger is currently used in three different 
contexts within the HETE project: 

• The Chicago ground location pipeline [7] uses 
spif fy-trigger to identify the burst sample 
time with maximal signal-to-noise in the WXM 
data. This sample is used throughout the subsequent 
location analysis. 

• A robot script that runs after every downlink 
uses spif fy-trigger to search for untriggered 
bursts in FREGATE band C (40-300 keV) 1.3s res- 
olution survey data. During normal HETE opera- 
tion, it tends to see about 1 possible GRB per week, 
above and beyond detecting all triggers picked up 
in flight that are sufficiently hard, and long (or short 
but bright) to register at this timescale and in this 
energy band. Figure 2 shows an example of such an 
event, which was confirmed by BeppoSax. 

• The general untriggered burst search described by 
Butler & Doty [8] uses spif fy-trigger in par- 
allel to Butler & Doty's wavelet trigger, and runs on 
all survey data products. GRB011212 was in fact 
detected on the ground in this pipeline, by both the 
wavelet algorithm and by spif fy-trigger. 



CONCLUSIONS 

The spif fy-trigger algorithm can probe a wide 
spectrum of burst timescales. It is still possible that ini- 
tialization with a very short might miss a very long, 
slow-rising event, or that a very long initial T^ u might 
cause the SNR of a weak, short event to be too diluted to 
register before convergence is reached. However, careful 
choice of the range of t^ u spanned by the initial simplex 
can address this issue to a large extent. 

In any event, the algorithm may be re-run with rad- 
ically different initial values of T^ u . For example, re- 
running the algorithm three times, with set initially 
to 0.1s, 3s, and 100s — with suitably chosen initial sim- 
plices — one may probe a range of timescales that would 
probably require hundreds of criteria for a traditional 
trigger algorithm to examine. 

In principle, there is no reason the 
spif fy-trigger algorithm could not be deployed in 
flight in a future mission. The floating-point operations 
that it performs are not particularly expensive, particu- 
larly for modern space computing hardware. While more 
complex than a traditional trigger, it is not vastly more 
so, and its complexity is offset by its great flexibility, 
configurability, and dynamic range of burst timescales to 
which it is sensitive. 
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